A form of the conservation equations for fluid dynamics is presented, deducedusing slightly less restrictive hypothesis than those necessary to obtain theWestervelt equation. This formulation accounts for full wave diffraction,nonlinearity, and thermoviscous dissipative effects. A two-dimensional finitevolume method using the Roe linearization was implemented to obtain numericallythe solution of the proposed equations. In order to validate the code, twodifferent tests have been performed: one against a special Taylor shock-likeanalytic solution, the other against published results on a High IntensityFocused Ultrasound (HIFU) system, both with satisfactory results. The code,available under an open source license, is written for parallel execution on aGraphics Processing Unit (GPU), thus improving performance by a factor of over60 when compared to the standard serial execution finite volume code CLAWPACK4.6.1, which has been used as reference for the implementation logic as well.
展开▼
